Bypassing Cowling's theorem in axisymmetric fluid dynamos 
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We present a numerical study of the magnetic field generated by an axisymmetrically forced 
flow in a spherical domain. At small enough Reynolds number, Re, the flow is axisymmetric and 
generates an equatorial dipole above a critical magnetic Reynolds number Rm c . The magnetic field 
thus breaks axisymmetry, in agreement with Cowling's theorem. This structure of the magnetic 
field is however replaced by a dominant axial dipole when Re is larger and allows non axisymmetric 
fluctuations in the flow. We show here that even in the absence of such fluctuations, an axial dipole 
can also be generated, at low Re, through a secondary bifurcation, when Rm is increased above 
the dynamo threshold. The system therefore always find a way to bypass the constraint imposed 
by Cowling's theorem. We understand the dynamical behaviors that result from the interaction of 
equatorial and axial dipolar modes using simple model equations for their amplitudes derived from 
symmetry arguments. 

PACS numbers: 47.65.-d, 52.65.Kj, 91.25.Cw 



It is strongly believed that magnetic fields of planets 
and stars are generated by dynamo action, i.e., self gen- 
eration of a magnetic field by the flow of an electrically 
conducting fluid Planets and stars being rapidly ro- 
tating, axisymmetric flows about the axis of rotation have 
often been considered in order to work out simple dy- 
namo models [2]. A major setback of the subject fol- 
lowed the discovery of Cowling's theorem, which stated 
that a purely magnetic field cannot be maintained by dy- 
namo action [3]. However, it has been shown that mag- 
netic fields with a dominant axisymmetric mean part can 
be generated when non- axisymmetric helical fluctuations 
are superimposed to a mean axisymmetric flow This 
has been recently observed in the VKS experiment A 
strongly turbulent swirling von Karman flow driven by 
two counter-rotating coaxial disks in a cylindrical con- 
tainer self-generated a magnetic field with a dipole mean 
component along the axis of rotation. This has been as- 
cribed to an alpha effect due to the helical nature of the 
radially ejected flow along the two impellers @. In this 
letter, we show that there exists another mechanism for 
bypassing the constraint imposed by Cowling's theorem, 
without the help of non axisymmetric turbulent fluctua- 
tions. The mechanism is as follows: the primary dynamo 
bifurcation breaks axisymmetry in agreement with Cowl- 
ing's theorem. Then, the Lorentz force generates a non 
axisymmetric flow component which can drive an axisym- 
metric magnetic field through a secondary bifuraction. 
We show that direct numerical simulations confirm this 
scenario and that the two successive bifurcation thresh- 
olds can be very close in some flow configurations. The 
existence of two competing instability modes, the axial 
and equatorial dipoles, can lead to complex dynamical 
behaviors. Using symmetry arguments, we write equa- 
tions for the amplitude of these modes that are coupled 



through the non axisymmetric velocity component. We 
show that the observed bifurcation structure and the re- 
sulting dynamics can be understood in the framework of 
this simple model. 

We first numerically integrate the MHD equations in a 
spherical geometry for the solenoidal velocity v and mag- 
netic B fields, 
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In the above equations, p is the density, p is the mag- 
netic permeability and a is the conductivity of the fluid. 
The forcing is f = fo F, where = s 2 sin(7r sb) , F z = 
e sin(7rsc), for z > 0, using polar coordinates (s, 0, z) 
(normalized by the radius of the sphere a) and oppo- 
site for z < 0. F^ generates counter-rotating flows in 
each hemisphere, while F z enforces a strong poloidal 
circulation. The forcing is only applied in the region 
0.25a <| z |< 0.65a, s < so. In the simulations pre- 
sented here, so = 0.4, b~ x = 2s$ and c _1 = s$. This 
forcing has previously been introduced to model the me- 
chanical forcing due to co-axial rotating impellers used 
in the Madison experiment [7[. Although performed in 
a spherical geometry, this experiment involves a mean 
flow with a similar topology to that of the VKS exper- 
iment. Such flows correspond to 52 + £2 flows in the 
Dudley and James classification [2j], i.e. two poloidal ed- 
dies with inward flow in the mid-plane, together with two 
counter-rotating toroidal eddies. We solve the above sys- 
tem of equations using the Parody numerical code |9j]. 
This code was originally developped in the context of the 
geodynamo (spherical shell) and we have here modified 
the code to make it suitable for a full sphere. We use 
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FIG. 1: Comparison of the magnetic field generated by an 
axisymmetric s2 + £2 flow in different geometries. Isovalue of 
the magnetic energy in a cylinder (from H) or a sphere. 




Reynolds number Re 
FIG. 2: Stability curves Rm c — f(Re) obtained with direct 
numerical simulations. In red: onset of the m — 1 (equa- 
torial dipole) dynamo mode; black: non-linear threshold of 
the m — (axisymmetric) mode; in blue the turbulent mode 
emerging from velocity fluctuations. The corresponding mag- 
netic structures are represented using magnetic field lines. 

the same dimensionless numbers as in [7], the magnetic 
Reynolds number Rm = /iocramax(|v|), and the mag- 
netic Prandtl number Pm = v\±§a. The kinetic Reynolds 
number is then Re = Rm/ Pm. 

The dynamo threshold Rm c is displayed as a function 
of Re in Fig. O Negative Re corresponds to a flow that 
is reversed compared to the VKS configuration, i.e. di- 
rected from the impellers to the center of the flow volume 
along the axis and radially outward in the mid-plane. For 
small enough i?e, the flow is laminar and axisymmetric. 
A magnetic field with a dominant equatorial dipole mode 
m = 1 is generated first (red curve in Fig. [2]). The ge- 
ometry of the field is displayed in the left inset of Fig. [2] 
and breaks axisymmetry as expected from Cowling's the- 
orem. This dynamo mode is similar to that obtained in 
cylindrical geometry, as illustrated in Fig [TJ 

For Re larger than about 300, the flow becomes tur- 
bulent and the equatorial dipole is then replaced by a 
dominant axisymmetric mode m = 0. Its threshold in- 
creases with Re in the parameter range of the simulations 
(blue curve in Fig. [2] and its geometry is shown in the 
right inset). These results are in agreement with [7(. It 
is remarkable that the axial dipole observed in the VKS 
experiment and ascribed to non axisymmetric fluctua- 
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FIG. 3: Bifurcation diagram of different magnetic modes 
varying Rm with fixed Re — +122. Error bars indicate the 
amplitude of oscillations. All other magnetic modes are very 
small compared to these ones. 

tions [6] can also be obtained in the present simulations 
even though the level of fluctuations is much smaller (the 
parameter range realized in the experiment being, by far, 
out of reach of present computer models). 

In addition, an axisymmetric magnetic field can also 
be generated at very low Re through a secondary bifur- 
cation from the equatorial dipole when Rm is increased. 
This corresponds to the black curves in Fig.[2j The corre- 
sponding mode is shown in the top left inset of Fig.[2j Bi- 
furcation diagram of Fig.[3]helps to understand the mech- 
anism by which this axisymmetric magnetic field is gen- 
erated. One can observe that the equatorial dipole first 
bifurcates supercritically for Rm = 88 when Re = +122. 
The back-reaction of the Lorentz force is twofold. First, it 
inhibits the axisymmetric velocity field, which decreases 
(orange curve in Fig. [3]). Second, and more importantly, 
it drives a non- axisymmetric m = 2 velocity mode (blue 
curve in Fig. [3]). Once the intensity of this flow be- 
comes strong enough, it yields a secondary bifurcation 
of the axisymmetric m = field mode. This is achieved 
for Rm = 205 (black curve in Fig. [3J. The ampli- 
tude of the equatorial dipole decreases immediately after 
this secondary bifurcation. We observe that the m = 
mode vanishes at higher Rm and then grows again above 
Rm = 425. Although the amplitude of the equatorial 
and axial modes behave in a complex manner as Rm is 
increased, we observe that they are anti-correlated, thus 
showing that they inhibit each other through the non- 
linear couplings. 

For Re < 0, Fig. [2] shows that the primary and sec- 
ondary bifurcations occur in a much narrower range of 
R m . The equatorial dipole mode is then close to marginal 
stability when the axial one bifurcates, and their non- 
linear interactions leads to complex time dependent dy- 
namics close to threshold as displayed in Fig. [U The 
equatorial mode (red curve) is generated first and satu- 
rates, but it drives the axial mode (black curve) through 
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FIG. 4: Time recordings of the total energy of the equatorial 
and axial dipolar modes, and of the m — 2 velocity mode for 
Re -76 and Rm 170. 
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FIG. 5: Bifurcation diagram of different magnetic modes 
varying Rm at fixed Re — — 76. Error bars indicate the am- 
plitude of oscillations. Dashed lines correspond to the maxi- 
mum values of the different modes in the relaxation regime. 
All other magnetic modes are very weak compared to these 
ones. 



the non axisymmetric part of the velocity field. The ax- 
ial dipole then inhibits the equatorial one that decays 
almost to zero. As a result, the flow is no longer driven 
away from axisymmetry by the Lorentz force. The ax- 
ial dipole thus decays and the process repeats roughly 
periodically. We observe that during one part of the cy- 
cle, the magnetic field is almost axisymmetric. It in- 
volves a strong azimuthal field together with a large ver- 
tical component near the axis of rotation, i.e. an axial 
dipole (see the left inset of Fig. [2]). These relaxation os- 
cillations, present only in the Re < case, occur only 
slightly above the threshold of the secondary bifurcation 
of the m = mode. Their period first decreases when 
Rm is increased, but then increases showing a divergence 
when the relaxation oscillations bifurcate to a stationary 
regime, as displayed in figure Above this transition, 
we observe bistability with the coexistence of two solu- 
tions: a nearly equatorial dipole, with a strong equatorial 



component and a weak axial one (labeled Mi in figure [5]) 
and a nearly axial dipole (labeled M2). 

We will show next that this competition between equa- 
torial and axial modes, and the resulting dynamics, can 
be understood using a simple model for the amplitudes 
of the relevant modes. We thus write 

B(r, t) = A(t) D eg (r) + c.c. + B(t) D asc (r) + • • • , (3) 

where D eg (r) (respectively D a£C (r)) is the eigenmode re- 
lated to the equatorial (respectively axial) dipole. A is 
a complex amplitude, its phase describes the angle of 
the dipole in the equatorial plane and c.c. stands for the 
complex conjugate of the previous expression. B is a real 
amplitude. As said above, the equatorial dipole (m = 1) 
generates a non axisymmetric flow through the action of 
the Lorentz force. The later depends quadratically on 
the magnetic field, this non axisymmetric velocity mode 
of complex amplitude V(t) thus corresponds to m = 2. 
Using symmetry arguments, i.e., rotational invariance 
about the z-axis which implies the invariance of the am- 
plitude equations under A — > Aex.pix,V — > Fexp2ix, 
and the B — > — B symmetry, we get up to the third order 

A = fiA- VA-a 1 \A\ 2 A-a 2 \V\ 2 A-a 3 B 2 A, (4) 
V = -vV + A 2 - /?i|A|V - f3 2 \V\ 2 V ~ PsB 2 V ,(5) 
B = -A£- 7l |A| 2 £ + 7 2 |F| 2 £-7 3 £ 3 . (6) 

\i is proportional to the distance to the dynamo thresh- 
old. Clearly v > 0, since the flow is axisymmetric below 
threshold. The coefficients of the quadratic terms can be 
scaled by an appropriate choice of the amplitudes. The 
term A 2 represents the forcing of the non axisymmetric 
flow by the Lorentz force related to the equatorial dipole. 
V A means that rotational invariance for the equatorial 
dipole is broken as soon as a non axisymmetric flow is 
generated. We have fixed its sign so that the bifurcation 
of the equatorial dipole remains supercritical Vai > 0. 
The equations for A and V (with B = 0) are the nor- 
mal form of a 1 : 2 resonance [10[ and have been studied 
in details in other contexts. In particular, it is known 
that this system can undergo a secondary bifurcation for 
which the phase of A begins to drift at constant veloc- 
ity when \i reaches a value such that | A \ 2 = 2 1 V | 2 . This 
corresponds here to a rotating dipole, at constant rate, in 
the equatorial plane. Consider now the equation for the 
amplitude B of the axial magnetic field. Taking A > 
and 73 > ensures that it cannot be generated alone, 
in agreement with Cowling's theorem. The term | V | 2 ^ 
describes the possible amplification of B from the non 
axisymmetric velocity field provided that 72 > 0. Al- 
though the system of amplitude equation dU|5]|6|) cannot 
be derived asymptotically from (jHEJ), it reproduces the 
phenomenology observed with the direct simulations for 
both signs of Re: when \i is increased, we either obtain 
relaxation oscillations as for Re < (parameters of fig. [6]) 
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FIG. 6: Numerical integration of the amplitude equations J2J- 
[5]|6]). Left: time recordings of the amplitudes of the equa- 
torial and axial magnetic modes interacting through the non 
axisymmetric velocity mode ( ji = l,ai = 0.3,0^ = 4, v — 
0.5, (3 2 = 1.5, A = 1.8, 7i = 0.5, 72 = 10, 73 = 0.5, all other co- 
efficients being zero.) Right: relaxation cycle in phase space 
involving the three unstable fixed points (real part of A, V). 



or a secondary bifurcation of the axial field as for Re > 
(same parameters with 72 = 1). The relaxation oscilla- 
tions are displayed in Fig. [6] (left). The model helps to 
understand the qualitative features observed in the di- 
rect simulation: it involves a solution corresponding to an 
equatorial dipole (Aq, Vb, B = 0) that can bifurcate to a 
mixed mode (Ai, Vi,Bi) involving a non zero axial field. 
In addition, two types of mixed modes can exist, one with 
a dominant equatorial dipole, say M\ = (Al, Vi, B\), and 
another with a dominant axial dipole M2 = (A2, V2, B2). 
Depending on the stability of these two solutions, we ob- 
serve either one of the mixed mode (depending on initial 
conditions), or a relaxation oscillation slowing down in 
the vicinity of these unstable fixed points and the origin. 
The system thus has three fixed points with both stable 
and unstable directions: the origin where both modes are 
zero, a point with a dominant equatorial dipole and a 
point with a dominant axial dipole. This situation leads 
to a heteroclinic cycle connecting these three unstable 
equilibrium points and corresponds to the relaxation os- 
cillations (see fig. [6j right) . 

In sodium flows driven by an axisymmetric forcing, 
such as the ones used in the VKS [5[, Madison and Mary- 
land experiments [11], one expects a possible competition 
between equatorial and axial dynamo modes. Indeed, 
the mean flow, if it were acting alone, would generate an 
equatorial dipole in agreement with Cowling's theorem. 
Our direct simulations show that a fairly small amount 
of non axisymmetric fluctuations (compared to the ex- 
periments) is enough to drive an axial (m = 0) dipole 
as observed in the VKS experiment for the mean mag- 
netic field. In addition, we show here that even without 
turbulent fluctuations, the non axisymmetric flow driven 
by the Lorentz force related to the equatorial dipole, can 
generate the axial one through a secondary bifurcation. 
The equatorial dipole can easily rotate in the equatorial 



plane, thus averaging to zero. The axial dipole then be- 
comes the dominant part of the mean magnetic field. 

It is striking that this mechanism that generates an 
axial dipole occurs much closer to the dynamo threshold 
when we go from the Re > to the Re < flow con- 
figuration, thus when the product of the helicity times 
the differential rotation is changed to its opposite value. 
For Re < 0, the shear layer in the mid-plane becomes 
favorable to an a — uo dynamo as soon as the axisymme- 
try of the flow is broken. For Re > 0, the flow near the 
impellers can play a similar role but the effect is weaker. 
This opens interesting perspectives for flows that can be 
used for future dynamo experiments: an a — uo effect 
driven by the strong vortices present in the shear layer 
close to the mid-plane can be favored by the Re < 
configuration. To wit, one can use either the optimized 
set-up described in [6| or propellors with the appropriate 
pitch in the VKS or Madison experiments. 

A competition between equatorial and axial dipolar 
modes could also account for secular variations of the 
Earth magnetic field. It would be interesting to check 
whether some features can be described with a low di- 
mensional model similar to the one used in this study. 

Computations were performed at CEMAG and IDRIS. 
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